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Abstract 

We calculate the gluon structure function of a color dipole in a new approach evaluating the matrix elements of SU(2) 
' 'gluon field operators separated along a direction close to the light cone. As vacuum state in the pure glue sector, we use 
t-H a variational ground state of the near-light-cone Hamiltonian. With a mean momentum fraction of the gluons fixed to 
the " experimental value" in a proton, the resulting gluon structure function for a dipole state with four links is compared 
^ qualitatively to the NLO MRST 2002 parameterization at Q 2 = 1.5 GeV 2 . 
O ■ 
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' — "i In Euclidean lattice simulations, using the operator 
'■^ product expansion the lowest moments of the meson and 
i" 1 nucleon structure functions have been evaluated [Il-Q- On 
Qnthe other side, Loop-loop correlation functions of tilted 
^ Wegner- Wilson loops computed in Euclidean space 0, Q 
t" 1 1 ban be related to the gluon distribution function of a color 
dipole which consists of a static quark and antiquark pair 
'connected by a Schwinger string. Alternatively, an ap- 
proach based on light cone dynamics @, H| refers to con- 
^vq tetituents moving along the light cone, as suggested by the 
f^i picture of Feynman scaling (as zeroth approximation). In 
t-H the light cone approach the non-perturbative QCD vac- 
__, .uum structure is hard to achieve within the Fock repre- 
sentation of free fields acting on a trivial vacuum. As we 
ON will demonstrate, the confining nontrivial QCD vacuum 
is essential to generate the correct interaction of colored 
constituents moving along light like trajectories. 

Therefore we have developed a ncar-light-cone (nlc) ap- 
^ proach in which we exploit the lattice formulation given 
C$ for SU (2) gluodynamics, benefitting from simplifications 
emerging in the light cone limit. In Ref. |9| we have con- 
structed a ground state wave functional of the nlc Hamilto- 
nian which is simpler than the ground state in equal-time 
theory. In Ref. [10| we have outlined the formalism to de- 
termine the gluon distribution function of a color dipole 
with this ground state wave functional. The present letter 
gives the main new results which one obtains following this 
approach. 

The gluon structure function g{xs) is the probability 
that a gluon carries a fraction of the longitudinal mo- 
mentum of the fast moving hadronic target. In light cone 
coordinates, it is given by the Fourier transform of the ma- 
trix element of the two-point operator Gi c (z~,z± ; 0, z±) 
of longitudinally separated gluon field strength operators 



X 



in a hadron state \h(p_, OjJ): 

1 1 f°° - 1 

9{xb) = — ^ dz~d 2 z ± e- 1XB P- z — 
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h(p-,6±) Gi c (z ,z±;0,z±) /i(p_,0jj 



with 



G lc (z-J ± ; 0,zl) =Y,Fl k {z~ J ± ) 

k=l 

x S^(z~, z ± ; 0, z ± ) F» k (0, z ± ) . (2) 



The hadron |/i(p_,0jj) is centered in transversal config- 
uration space at b± = 0± and carries a longitudinal mo- 
mentum p_ . The index "c" indicates that the connected 
matrix element has to be taken. The Schwinger line 
SA(z~,z± ; 0, z±) in the adjoint representation and run- 
ning along a light like path is inserted between the gluon 
field strength operators F^_ k (Q,z±) and F^ k (z~ , z±). The 
importance of the Schwinger lines along the light cone has 
been demonstrated e.g. in the loop-loop correlation model 
where hadron-hadron scattering cross-sections are calcu- 
lated from Wegner- Wilson loop correlation functions @. 
In another language, the eikonal phases arising from the 
strings along the x~- direction describe "final state" inter- 
action effects which distinguish structure functions from 
parton probabilities 

We are using near-light-cone coordinates which allow 
us to implement light front quantization as a limit of equal 
time quantization. The definition of the temporal nlc 
coordinate x + contains an additional external parameter 
rj which facilitates a smooth interpolation between equal 
time quantization (rj = 1 , x + = x ) and light cone quan- 



tization (rj = 



r+ - 



= 1/2 (x° + x 3 )). The definitions 



Preprint submitted to Elsevier 



February 4, 2010 



arc: 



1 



T 



'1 



r, 2 )x 3 



The r\ — > limit can be interpreted as the infinite momen- 
tum frame limit in which the partons of the color dipole 
move with infinite momentum. The use of nlc coordinates 
has the advantage that no quantum constraint equations 
have to be solved. 

In the SU (2) lattice formulation we have two transver- 
sal gauge fields Aj (j = 1, 2) and one longitudinal one A_ 
(J = — ) that are represented by the corresponding link 
matrices Uj{x). The A + component of the gauge field is 
set equal to zero. As a result one has the Gauss-law con- 
straining the entire Hilbert space to the physical sector of 
gauge invariant states. The gluon dynamics is determined 
by the Hamiltonian which has been derived in Ref. It 
represents the gluon energy density on the lattice. The 
QCD coupling constant enters as A = 4/g 4 in the SU(2) 
case, 

H — ^l47s?0? nl( ' f)2 + 




Figure 1: The generalized Wegner- Wilson loop generating the matrix 
element between the color dipole states. Quark q and antiquark q 
are connected by n transversal links. Only one transversal direction 
is shown. The insertion represents the gluonic two-point operator 
with an electric and a magnetic field operator. 



operator in Eq. ([2]) has to be replaced by 



i 2 

G lc (z-,z ± ; 0,z x ) = -J2 {F* k (z-,z ± ) 



xS^ b (z-,z x] 0,z x )U b k (0,z x ) + {Fl k ^U%}) . (5) 

In Ref. [9] we have determined a variational gluonic 
ground state wave functional |\&n) which consists of a prod- 
uct of single-plaquette wave functionals with two variation- 
ally optimized parameters po — po(A, rj) and So = <5o(A, 77), 

|*o) = *o[U] |0> = y/Ifye'W |0) , 

x lfe=l J 

th^* * s a normanza ti° n factor. Here, the state |0) represents 
the trivial ground state which is annihilated by the field 
momenta 11^ (2?) conjugate to the links, 

II£(z) |0) = and (0| H%(x) = for Vf,fc,a. (6) 

We have optimized this ansatz and extrapolated the pa- 
rameters po^o to the light cone 77 — ^ (see Ref. Q). 

Since the nlc Hamiltonian in Eq. §3§ contains only 
gluon fields, we cannot derive a full hadronic wave function 
from this Hamiltonian. We have to make a model taking 
care of the gluon structure alone while treating the quarks 
schematically. Our model represents a dipole localized in 
transversal configuration space at a fixed center of mass 
position b± = 0. The distance d± between quark and an- 
tiquark is bridged by a Schwinger line along some path 
C± in the transversal plane. The in and out dipole states 
form a Wegner- Wilson loop in the eikonal approximation 
as sketched in Fig. [TJ 

The longitudinal lattice momenta must be n-fold (in- 
teger) multiples of 2ir/N-, with < n < iV_/2 — 1, since 
the longitudinal light cone momentum for an on shell par- 
ticle is always positive. The momentum p_ of the target is 
chosen (in lattice units) as the largest momentum in order 
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Here, U12 are purely transversal plaquettes, and U-k are 
longitudinal-transversal plaquettes. 

The Hamiltonian contains the chromo-electric field strcnj 
operators Uf(x), which are canonically conjugate to the 
links Ui{x). They obey commutation relations which fol- 
low from the corresponding continuum relations, 



[IlUx),U 3 m - (a a /2) r,i.f)6,,, ; S, 



(4) 



The constant f\ in Eq. is the product fj — rj ■ £ of the 
near-light-cone parameter r\ and an eventual anisotropy 
parameter £ = a_/ai, the ratio of lattice spacings in lon- 
gitudinal and transverse directions. If one chooses 77 = 1 
and varies £, one simulates an anisotropic equal time the- 
ory. In the limit £ — ¥ one ends up with a system, 
which is contracted in the longitudinal direction. Ver- 
linde and Verlinde [l2| and Arefeva [l3| have advocated 
such a lattice to describe high energy scattering. A lon- 
gitudinally contracted system means that even the mini- 
mal momenta become high in longitudinal direction. The 
limit £ — > leads to the same physics as the light cone 
limit rj — > with isotropic lattice spacings in longitudinal 
and transverse directions. In both cases the nlc Hamil- 
tonian is dominated by the terms proportional to l/ij 2 
involving transverse chromo-electric and chromo-magnetic 
fields. For a nlc-Hamiltonian formulation, the two-point 
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to allow for the maximum resolution in the gluon distribu- 
tion function [14| . The longitudinal lattice gluon momenta 
then have the resolution 

Ap 9 _ 2 



P- = 



2ir ,N. 



N- 



P- 
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(7) 



In order to have a high resolution, the extension of the 
lattice in the longitudinal direction has to be very large. 

We impose the quark dynamics of the color dipole ex- 
ternally. Since the total hadron longitudinal momentum 
is given by the sum of the momenta of its constituents 
( quarks + gluons), the typical mean gluon momentum is 
taken from experiment. For a rough qualitative compar- 
ison, we use the Mi?5T-parameterization [l5| of the pro- 
ton SU{2>) gluon distribution function at the input scale 
Q 2 w 7r 2 /a^ = 1.5 GeV 2 corresponding to A w 10, and 
assign a mean momentum fraction p^_ — 0.38 p_ to the 
gluon system of the SU(2) color dipole. 

The gluon distribution function for a one-link dipole 
P- 9i (p 9 - ; P- ) with total gluonic momentum p_ has been 
computed on lattices with iV_ = 20, 30, 50, and 100 sites 
in the longitudinal direction. Fig. [2] demonstrates the ef- 
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Figure 2: Gluon distribution function p 9 _ gi(p B _ 



0.38 p_) for 



a single transversal link, n = 1. The dependence on the longitudinal 
lattice size 7V_ is shown for a lattice coupling A = 4/g 4 = 10 (see 



Eq. (4)). The average gluon momentum 
adjusted to (xg) p_ = 0.38 p— (see text). 



p_ was input and 



feet of increasing the number of longitudinal lattice sites, 
i.e. approaching the infinite volume limit. Scaling for 
A = 10 seems to be obeyed for longitudinal lattice ex- 
tensions larger than iV_ = 50. Realistic lattice simulations 
need quite large longitudinal lattice sizes. The smearing of 
the distribution function is due to the gluon dynamics in- 
corporated in the Wegner- Wilson loop expectation value. 
Thus, the area law behavior of the Wegner- Wilson loop 
yields a non-trivial gluon wave function which broadens 
the distribution. 

For larger A, i.e. smaller QCD gauge coupling g 2 , the 
peak in the one-link distribution function becomes nar- 
rower. In the extreme weak coupling limit, when the link 



reduces to a single gluon, the gluon distribution function 
becomes totally sharp. On the other hand, for strong 
coupling one has a broad momentum distribution peaked 
around p_ . 

The one-link dipole gluon distribution is the basic build- 
ing block from which the multiple link dipole gluon distri- 
bution function of a hadron can be constructed. The ac- 
tual hadronic state arises from a superposition of multiple 
link configurations. Let the wiggly strings S q q (c.f. Fig.QJ 
connecting the quark and antiquark have a fixed number 
of transversal links since all the dipole configurations with 
fixed transversal length n have the same energy. Then one 
rotates the hadron in the transversal plane by summing 
over randomly chosen curves C±, in order to project this 
state on angular momentum J z — 0. ^From the random 
walk follows that for n-links , the hadron has an average 
radius squared R\ proportional to n: Hence, the area of 
the hadron scales with the number of links: 
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(8) 



Due to small Sq in the light cone limit, the ground state 
wave functional allows for ordinary strong coupling meth- 
ods. This implies especially that one needs incoming and 
outgoing states sharing the same transversal links connect- 
ing the quark and antiquark. Therefore, the gluon distri- 
bution can be obtained for a special string elongated along 
only one of the transversal axes (c.f. Fig. [T]). 

The computation of the n-link gluon distribution func- 
tion g n is done in analogy to the computation of the sin- 
gle link gluon distribution. In strong coupling the total 
loop factorizes, therefore the n-link distribution function 
is given by the product of a splitting function P„_>. n _i mul- 
tiplying the gluon distribution function with n — 1 links. 
In the emerging recursion relation all possible intermediate 
momenta of the substring are summed over: 



(p g -; P s -) 



2tt x ^ 

N- ^ 



9n-l 



(p 9 -;p s -') 



Pn- 



l(p*,p£') 



(9) 

The splitting function P n ^ n -i{P-,P-') denotes the prob- 
ability that a string with n transversal links and total mo- 
mentum p^_ splits into a string with n—1 transversal links 
and total momentum p^_' and is given in Ref. (lp| . The 
initial condition for the recursion relation Eq. ^ is given 
by the one- link dipole function gi(p 9 _',p^_). The computa- 
tion is purely arithmetic for small So, and we can use a 
large longitudinal lattice with N- = 1000 lattice sites. If 
one increases the number of transversal links, the gluons 
have access to a larger region in phase space due to the 
splitting function P„^„_i in Eq. ©. Therefore the total 
gluon momentum will be partitioned among more gluons. 
Hence, it becomes more likely to find a gluon with a small 
fraction of the total momentum. 

In Fig.[3l we show the theoretical gluon structure func- 
tion for a n = 4 SU(2) link dipole as a function of the gluon 
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Figure 3: Qualitative comparison of the gluon distribution function 
%B9n{zB',P—) for n = 4 transversal SU(2) links with the MUST 
gluon distribution function of a proton at Q 2 = 1.5 GeV 2 . The av- 
erage gluon momentum was taken as input for the lattice calculation 
on a JV_ = 1000 lattice at gauge coupling A = 10. 
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fractional momentum xb = P-/P-- To have a rough qual- 
itative comparison, we also show the MRST gluon distri- 
bution function at Q 2 = 1.5 GeV 2 which is for real pro- 
tons and SU(3). As before, the first moment of the lattice 
gluon distribution function has been fixed in this figure to 
the value (x B ) = 0.38 at Q 2 = 1.5 GeV 2 . We choose four 
links to be consistent with the size of a proton and the re- 
lation (Rj_) = na\/2 and a transversal lattice spacing of 
a± R3 0.65 fm. The lattice gluon distribution function for 
a color dipole shows a similar behavior as the phenomeno- 
logical MRST-gluon distribution function for a proton. It 
has been shown in the loop-loop correlation model [lit 
that it is a quite good approximation to consider the three 
valence quarks in the proton to be arranged as a dipole 
with a diquark and a quark. 

The model presented here also shows that xb 9{xb) 
for the gluon at small xb becomes proportional to the 
hadronic size R\. This coincides with the empirical soft 
Pomeron behavior of hadronic cross sections. To evolve 
the structure function with increasing resolution Q 2 and 
with decreasing xb in Hamiltonian lattice QCD, one needs 
a more sophisticated ground state wave functional respect- 
ing scaling with the lattice spacing. 
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